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Abstract. We compute theoretical nova rates as well as type la SN rates in galaxies of different morphological 
type (Milky Way, ellipticals and irregulars) by means of detailed chemical evolution models, and compare them 
with the most recent observations. The main difference among the different galaxies is the assumed history of star 
formation. In particular, we predict that the nova rates in giant ellipticals such as M87 are ~ 100 — 300 nova yr~^, 
about a factor of ten larger than in our Galaxy (~ 25 nova yr~^), in agreement with very recent estimates from 
HST data. The best agreement with the observed rates is obtained if the recurrence time of novae in ellipticals is 
assumed to be longer than in the Milky Way. This result indicates that the star formation rate in ellipticals, and 
in particular in M87, must have been very efficient at early cosmic epochs. We predict a nova rate for the LMC 
of 1.7 nova yr~^, again in agreement with observations. We compute also the K- and B-band luminosities for 
ellipticals of different luminous mass and conclude that there is not a clear trend for the luminosity specific nova 
rate with luminosity among these galaxies. However, firm conclusions about ellipticals cannot be drawn because 
of possible observational biases in observing these objects. The comparison between the specific nova rates in the 
Milky Way and the LMC indicates a trend of increasing nova rate passing from the Galaxy towards late-type 
spirals and Magellanic irregulars. 
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1. Introduction 

^ ' Classical novae are binary systems, where one component 
Jv>( , is a white dwarf (either carbon-oxygen or neon-oxygen- 

5^ magnesium), accreting hydrogen rich material from a less 
1 evolved companion. When the pressure and temperature 
at the bottom of the accreted layer exceed critical val- 
ues (e.g. Livio 1992) the nuclear reactions are ignited and 
novae exhibit a sudden and rapid increase of their bright- 
ness. The thermonuclear runaways are responsible for the 
production of interesting chemical species, such as "'^^C, 
i^N, "O, 22Na, 26 Al, 22Ne and possibly ^Li (see Jose' 
& Hernanz 1998, Gehrz et al. 1998, Delia Valle et al. 
2002) and for ejection of the accreted envelope. Although 
the matter processed through nova explosions is a small 
fraction (a few 10^'^) of the total mass of the interstel- 
lar gas and dust in the Milky Way, the concentrations of 
these rare isotopes in the nova ejecta can be enhanced, 
relative to solar abundances, by factors lO^^"^, therefore 
making novae non-negligible contributors to the galactic 
nucleosynthesis (Romano et al. 1999, 2001; Romano & 
Matteucci 2003). Novae have been observed not only in 
our Galaxy but also in the Local Group (e.g. Hubble 1929, 
Arp 1956 and Rosino 1964) and beyond, in Virgo (Pritchet 



& van den Bergh 1987) in M81 (Shara, Sandage & Zurek 
1999) and in Fornax (Della Valle & Gilmozzi 2002). From 
these studies some estimates of nova rates have been de- 
rived and they are shown in Table 1, where the Hubble 
type of the galaxies is indicated in Col. 1, the nova rate 
in yr~^ in Col. 2, the K-band luminosity of the object in 
Col. 3 and in Col. 4 the luminosity specific nova rate in 
units of yr^^lO~^°L0.K {LSNRk)- These rates are taken 
from Delia Valle (2002) and Shatter et al. (2000). 

In this paper we use the chemical evolution mod- 
els elaborated by Chiappini et al. (2001), Matteucci et 
al. (1998) and Calura et al. (2003) for the Milky Way, 
Ellipticals and late Spirals, respectively, to compute the 
nova rate in 'disk' dominated and 'bulge' dominated sys- 
tems with the intent of checking whether the different ob- 
served nova rates can be reproduced under the assumption 
that the main difference between galactic morphological 
types should be ascribed to the different histories of star 
formation. Moreover, our study is aimed at deriving im- 
portant constraints on some key parameters of nova the- 
ory, which cannot be directly inferred from observations, 
such as the recurrence time between the outbursts, the 
fraction of binary systems which give rise to nova sys- 
tems, the delay time with which the first novae appear 
relative to the beginning of star formation. The use of such 
models is appropriate since they can reproduce the basic 
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Table 1. Extragalactic nova rates and luminosity-specific nova rates from Delia Valle (2002). The latest estimate for 
M87(M87*)is from Shara & Zurek (2002). With Virgo E we refer to other elhpticals in the Virgo clusters. 



Galaxy Hubble type novae (yr~^) LK(lO"L0,jf) LSNRk (yr~^ (1O^°L0,k)~^) 



M31 


Sb I-II 


29± 4 


18. 5± 2.4 


1 ccj. on 


M33 


Sc II-III 


4.6± 0.9 


1.17±0 .15 


3.93±0.93 


M51 


Sbc I-II 


18±7 


16.6±3.4 


1.09±0.47 


M81 


Sb I-II 


24±8 


13.3±2.7 


1.80±0.71 


M87 


EO 


91±34 


39.8±8.2 


2.30±0.99 


M87* 




200 to 300 






MlOl 


Sc I 


12.0±4.0 


12.4±2.6 


0.97±0.38 


NGC5128 


SO(Pec.) 


28±7 


7.94±1.10 


3.52±1.01 


LMC 


SBm III 


2.5±0.5 


0.46±0.06 


5.42±1.30 


SMC 


Im IV-V 


0.7±0.2 


0.14±0.02 


2.19±1.49 


Virgo 


E 


160±57 


49.7±10.5 


3.22±1.33 


NGC1316 


SO(Pec.) 


90 to 180 


? 


? 



parameters characterizing the respective parent galaxies 
such as the rate of type la SNe, B and K luminosities of 
the parent galaxies, the [a/Fe] ratios and global metallic- 
ities in the stellar popiilations and gas. Finally, we note 
that a fraction of type la SNe can originate from WDs in 
CV-like systems, after reaching the Chandrasekhar mass 
limit (e.g. Livio 2000). Whether these systems can give rise 
to novae or supernovae depends crucially on the mass of 
the underlying white dwarf and mass accretion rate (see 
Nomoto et al., 1984). Therefore, when trying to predict 
nova rates in galaxies it is obvious to compute and try to 
reproduce the observed SNIa rates. The observed type la 
rates in galaxies of different morphological type are shown 
in Table 2. 

While the rates of novae are expressed in yr~^ 
or in LSNRk, the rates of SNIa are expressed 
in SNub (SN (100yr)-il0-i°LBQ) and SNuk, (SN 
(100yr)~^10~^°LKQ). The paper is organized as follows: 
in Sect. 2 we describe how to compute the nova and SN 
rates in galaxies and compare the results with observations 
and finally, in Sect. 3, we draw some conclusions. 



siders an almost independent evolution between the halo 
and thin-disk components (see also Pagel & Tautvaisiene 
1995), and it assumes a threshold in the star formation 
process, (see Kennicutt 1989, 1998; van der Hulst et al. 
1993; Martin & Kennicutt 2001). The model weU repro- 
duces the majority of observational constraints about the 
abundances of heavy elements both locally and in the 
whole disk. 

Table 2. SNIa rate (SNub) from the combined search 
sample of CappcUaro ct al. (1999) and SNIa rate in SNuk 
(Delia Vallc & Livio 1994) SNub and SNuk are in units 
ofSN (100yr)-ilO-i"LBQ and SN (100yr)-ilO-i°LKQ, re- 
spectively, for h = 0.75. 

Hubble type SNela {SNub) SNela {SNuk) 



E-SO 0.18±0.06 0.030 ±0.01 

Sa^Sb 0.18±0.07 0.045 ±0.015 

Sc-Sm 0.40±0.08 0.140 ±0.03 



2. Modelling the nova and supernova rates 

2.1. The Milky Way 

The model for the Galaxy assumes two main infall 
episodes for the formation of the halo-thick disk, and the 
thin-disk, respectively. The timescale for the formation of 
the thin-disk is much longer than that of the halo, im- 
plying that the infalling gas forming the thin-disk comes 
not only from the halo but mainly from the intergalactic 
medium. The timescale for the formation of the thin disk 
is assumed to be a function of the galactocentric distance, 
leading to an inside-out picture for the Galaxy disk build- 
up, according to the original suggestion of Matteucci & 
Frangois (1989). The two-infall model differs from other 
models in the literature mainly in two aspects: it con- 



The star formation rate (SFR) adopted here has the 
same formulation as in Chiappini et al. (1997): 

*<«.') = "<"(|Sy)""'(^)'""«:-(«.')« 

where v{t) is the efficiency of the star formation process, 
S(i?, t) is the total surface mass density at a given radius 
R and given time t, 'E{RQ,t) is the total surface mass 
density at the solar position, Ggas {R, t) is the surface gas 
density normalized to the present time total surface mass 
density in the disk E£)(i?, tcai), ^Gai = 13 Gyr is the age of 
the Galaxy, R© = 8 kpc is the solar galactocentric distance 
(see Reid 1993). 
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Galaxy Routburata{tGai) ~ 25yr \ as in Romano et al. 
(1999), for the following reasons. Given the fact that obser- 
vations of novae in nearby galaxies would avoid, or at least 
minimize, most of the problems encountered by Galactic 
observations, such as interstellar extinction in the Galactic 
disk. Shaftcr (1997) shows that the nova rate based on 
Galactic observations can be made consistent with the 
rate predicted from the extragalactic data (Delia Valle 
et al. 1994) particularly if the Galaxy has a strong bar 
oriented in the direction of the Sun. In this latter, most 
favourable case, the suggested value is near k. 20yr~^, oth- 
erwise, if the bar is weak or misaligned, the global rate can 
be w 30yr-i. 

Table 3. Results for the model for the Galaxy, IMF Scalo 
(1986), a = 0.01, ^=0.05, tGai=13 Gyr. 



The gas surface density exponent, k, is set equal to 
1.5, in order to ensure a good fit to the observational con- 
straints in the solar vicinity. This value is also in agree- 
ment with the observational results of Kennicutt (1998), 
and with N-body simulation results by Gerritsen & Icke 
(1997). The star formation efficiency is set to v = 2 Gyr"-*^ 
for the Galactic halo, whereas \i \s v = 1 Gyr~^ for the 
disk; this is to ensure the best fit to the observational 
features in the solar vicinity. The star formation rate be- 
comes zero when the gas surface density drops below a 
certain critical threshold. We adopt such a threshold den- 
sity to be {(Tgas)th ~ 4 Mq pc~^ for the Galactic halo, 
and {agas)th ~ 7 Mq pc~^ for the disk (see Chiappini 
et al. 2001). The assumption of such threshold densities 
naturally produces the existence of a hiatus in the SFR 
between the halo-thick disk phase and the thin-disk phase. 
This discontinuity in the SFR is observed in the [Fe/0] vs. 
[0/H] (Gratton et al. 2000) and in the [Fc/Mg] vs. [Mg/H] 
(Fuhrmann 1998) plots. The IMF is that of Scalo (1986), 
assumed to be constant in time and space. We compute 
the nova systems formation rate {Rnovae{i)) at the time 
t as a fraction of the formation rate of white dwarfs at a 
previous time t — At (see D'Antona & Matteucci, 1991): 

^.8.0 

Rnovae{t) = a tp{t ~ Tm - At)(j){m)dm. (2) 
Jo.8 

Here At is a delay time whose value has to be fixed to 
guarantee the cooling of the white dwarf (WD) to a level 
that ensures a strong enough nova outburst. We assume 
a distribution of At varying from 1 to 5 Gyr. It is worth 
noting that assuming a constant average delay of At=2 
Gyr, as assumed in previous papers (e.g. Romano et al. 
1999; 2001), produces negligible differences in the final 
results. 

It is assumed that all stars with initial masses in the 
range 0.8 - 8 M© end their lives as WDs. ip{t) is the SFR 
as defined in Eq. (1), is the lifetime of the star of mass 
m and (f>{ni) is the IMF. 

The rate of nova eruptions is related to the WD for- 
mation rate by: 

^outbursts {t) = (xRwDs {t)n, (3) 

where aRwDs = Rnovae{t) is the formation rate of WDs 
in binary systems, as defined in Eq. (2), which will give 
rise to nova eruptions and n = 10^ is the average number 
of nova outbursts for a typical nova system (Ford, 1978; 
Bath & Shaviv 1978; Shara et al. 1986). This number, 
which is just an average value, is inversely proportional 
to the recurrence time. The parameter a, which here is a 
constant in space and time, can be fixed by reproducing 
the rate of nova outbursts in our Galaxy at the present 
time. Unfortunately, estimates of this quantity in the cur- 
rent literature show a large spread. Predictions based on 
scalings from extragalactic nova surveys suggest low val- 
ues, whereas estimates based on extrapolations of Galactic 
nova observations give the highest rates. We consider for 
the present time observed rate of nova outbursts in the 



nova (yr-^) nova (LSNRk) SNela (100?/r~^) SNela (SNu) 



25 2.1 0.3 0.16 



Our first choice is to fix a = 0.01 (coupled with n = 

10^) so to reproduce the present observed nova rate of 
25yr~^ and to use the same a in the models for elliptical 
and irregular galaxies. This value of a is the same as in 
previous papers (e.g. Romano et al. 2001; Matteucci et al. 
1999). 

Finally, the SNela rate has been computed following 
Grcggio & Rcnzini (1983) (see also Matteucci & Greggio 
1986) and is expressed as: 

RsNela = A (j){MB) / /(M)V'(i " TM^rf/i dMB(4) 

J l^m 

where M2 is the mass of the secondary, Mb is the to- 
tal mass of the binary system, /i = M2/MB, jJ-m = 
max{M2{t)/MB,{MB - o!5Mb„)/Mb}, Mb^ = 3 M©, 
Mbm = 16 Mq. (f>{MB) is the initial mass function for 
the total mass of the binary system, /(/x) is the distri- 
bution function for the mass fraction of the secondary, 
/(/i) = 2^+^(1 + -f)ii'y, with 7 = 2; A = 0.05 is the frac- 
tion of systems with total mass in the appropriate range, 
which give rise to SNIa events. This quantity is fixed by 
reproducing the observed SNe la rate at the present epoch 
(see also Madau et al. 1998), and it provides excellent fits 
to the abundances and abundance ratios in the Galaxy. 

The results for the Galaxy are shown in Table 3 and 
Fig. 1. In Table 3 we report the predicted nova and SN la 
rate for the Milky Way, adopting an age of 13 Gyr, accord- 
ing to a cosmology with fi^ = 0.3, JIa = 0.7, h = 0.7 and 
Zf = 10. We computed the K- and B- band luminosity 
of the Milky Way by means of the spectro-photometric 
model of Jimenez et al. (1999): we find, for the present 
time, luminosities in very good agreement with the es- 
timated ones. The estimated Galactic blue luminosity is 



Prancesca Matteucci et al.: Modelling the nova rate in galaxies 



4 

Lbg,, ~ 2 • 10i°Lo (see van den Bergh 1988) and the K- 
luminosity is Lk^al ~ 12-10^°Lq (see Shatter el al. 2000). 
In Figure 1 wc show the behaviour of the star formation 
rate in the framework of the two-infall model, where the 
gap between the formation of halo-thick disk and thin-disk 
is evident, as well as the behaviour of the nova and SN la 
rates. These latter show only a small discontinuity in cor- 
respondance with the star formation gap and the reason 
is that they both depend upon the past star formation. 
The nova rate increases continuously up to the present 
time due to the long delay assumed for these systems (the 
lifetime of the WD plus the At), whereas the SNIa rate 
reaches a maximum and then it stays flat afterwards. 

2.2. Ellipticals 

In the models for ellipticals, the galaxies can be considered 
initially as an homogeneous sphere of gas with luminous 
masses in the range 10^ — 10^^ Mq. A single zone interstel- 
lar medium with instantaneous mixing of gas is assumed 
throughout. 

The adopted model parameters are shown in Tables 
4 and 5, where we show the initial luminous masses for 
the model galaxies, the assumed star formation efficiency 
parameter, the effective radius, the predicted time for the 
occurrence of a galactic wind taw and the predicted K- 
and B-luminositics. The parameter A is varied according 
to the assumed IMF in order to reproduce the present time 
observed type la SN rate (Matteucci & Gibson, 1995). On 
the other hand, wc kept the parameter a to be the same 
as in the Milky Way. The model is described in Matteucci 
et al. (1998) where we address the reader for more details. 

The star formation rate is given by: 

i>{t) = ypgas{t)/p{Q), (5) 

i.e. normalized to the initial total volume density. V'(0 is 
assumed to drop to at the onset of the galactic wind. 
The quantity v is expressed in units of Gyr^^ and repre- 
sents the efficiency of star formation, namely the inverse 
of the time scale of star formation. The star formation is 
assumed to stop after the development of a galactic wind 
occurring before than 1 Gyr, from the beginning of star 
formation, for all the galaxies listed above. Therefore, the 
star formation rate in these galaxies can be considered as 
a strong burst which does not last more than 1 Gyr and 
is shorter in more massive systems. This is obtained by 
assuming that the star formation efficiency increases with 
galactic mass thus producing an " inverse wind" effect, 
as described in Matteucci (1994), where the galactic wind 
occurs before in more massive than in less massive ellipti- 
cals. As a consequence of this, the star formation period 
is longer in smaller systems thus allowing the SNe la to 
substantially pollute the ISM. This effect can explain the 
observed increase of the [a/Fe] ratio with galactic mass 
(Worthey et al. 1992; Matteucci 1994), which is not ob- 
tained in a classic wind scenario, where the more mas- 
sive objects form stars for a longer time (Larson, 1974). 



The galactic wind develops as a consequence of the en- 
ergy transfer from SNe into the ISM. In fact, when the 
thermal energy of the gas becomes larger than the bind- 
ing energy of the gas, the wind starts (Arimoto & Yoshii, 
1987; Matteucci & Tornambe' 1987; Matteucci 1992, 1994; 
Pipino et al. 2002). In order to compute the binding en- 
ergy of the gas some assumptions have to be made about 
the galactic potential well. In particular, it is assumed that 
all ellipticals possess heavy but diffuse dark matter halos; 
a ratio between the half-light radius and the radius of the 
dark matter core Riuminous / Rdark=^AQ and a ratio dark 
to luminous mass of 10 are assumed. 

The B and K luminosities for the ellipticals of different 
masses are computed by means of the spectro-photometric 
model of Jimenez et al. (1999) and are used to compute 
the SN la rate in SNu and the nova rate per unit of Lk, 
respectively. 

The main differences between the model for the 
Galaxy, and the model for an elliptical galaxy concern 
the different SFR, which is much stronger in the earli- 
est stages, and then is set to zero after the galactic wind 
in the case of ellipticals. It is worth noting that galactic 
winds in the case of elliptical galaxies seem necessary to 
explain their lack of gas and the chemical enrichment of 
the ICM, whereas the Galactic model does not take into 
account the occurrence of a strong wind, mainly because 
of the strong gravitational potential well associated with 
our Galaxy, and the presence of the gas in the Galactic 
disk. 

We adopt two different initial mass functions (IMFs): 
the Salpeter {xjmf — 1.35) one and the Arimoto & Yoshii 
(1987) {xjMF — 0.95) one. In fact, successful models of 
chemical evolution of ellipticals have shown that the Scalo 
IMF is not suitable for these galaxies (e.g. Matteucci & 
Gibson, 1995). 

Table 6 shows the results of the model I (Salpeter IMF), 
whereas Table 7 shows the results of the model II (x=0.95 
IMF). These models have been computed by assuming the 
same a and the same n as for the Milky Way. 

In particular, in the first column we report the initial 
galactic luminous mass, in the second column the nova 
rate in units of yr^^, in Col. 3 the luminosity specific nova 
rate , in Gol. 4 the SN la rate in units of (lOOyr)^^ and 
in Col. 5 and 6 the SN la rate in units of SNub and 
SNuk, respectively. The predicted nova rates are quite 
large for massive ellipticals ranging between 300 and 1000 
nova yr~^ which is about a factor 3-10 larger than it is 
derived from observations. If these values were realistic, 
the discrepancy could be partially due to an observational 
bias affecting the ground based nova surveys due to poor 
spatial resolution and bright limiting magnitude. An indi- 
cation in this direction comes from the nova rate of M87 
recently provided by Shara & Zurek (2002) on the basis 
of HST archive images, which is a factor 2-3 larger than 
previous ground-based estimates (see Table 1). Another 
possibility is that in early type galaxies the recurrence 
time between two consecutive nova explosions is consid- 
erably longer than in late spirals as a consequence of the 
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Table 4. Model I: a; = 1.35, a = 0.01, A = 0.18, tgai = 13Gyr. 



M;„,„(M0) 


u (Gyr-i) 


R-e// 


(kpc) 


taw (Gyr) 




Lb (101°Lsq) 


10« 


2. 


0.5 




1.4 


0.038 


0.0091 


101° 


5. 


1.0 




1.3 


0.40 


0.093 


10" 


11. 


3.0 




0.6 


5.9 


0.70 


4-1011 


15. 


6.0 




0.4 


17.0 


3.9 


e-iQii 


16.5 


7.0 




0.4 


25.0 


5.9 


1012 


20. 


10.0 




0.4 


45.0 


9.6 


Table 5. Model II: x = 


0.95, a = 0.01 


, A = 0.05, tgai = ISGyr. 






I' (Gyr-i) 


Re// 


(kpc) 


taw (Gyr) 


Lk (101°Lkq) 


Lb (10i°Lsg) 


10^ 


2. 


0.5 




1.4 


0.034 


0.0044 


10i° 


5. 


1.0 




0.9 


0.42 


0.049 


10" 


11. 


3.0 




0.6 


4.3 


0.53 


4-10" 


15. 


6.0 




0.4 


19.0 


2.3 


6-10" 


16.5 


7.0 




0.4 


37.0 


3.2 


1012 


20. 


10.0 




0.3 


47.7 


6.2 



Table 6. Results for model I: nova and SN la rates. 

Mi„^(Mo) novae (yr'i) novae (LSNRk) SNe la (lOOyr-i) SNe la (SNus) SNe la {SNuk) 



10^ 0.6 15.8 0.0026 0.28 0.07 

IQi" 8 20.0 0.027 0.29 0.08 

10" 97 16.4 0.19 0.27 0.03 

4-10" 409 24.0 0.54 0.14 0.03 

6-10" 635 25.4 0.78 0.13 0.03 

10i2 1121 24.9 1.00 0.10 0.02 

Table 7. Results for model II: nova and SN la rates. 

Mi„^(M0) novae (yr'i) novae (LSNRk) SNe la (lOOyr'i) SNe la (SNub) SNe la {SNuk) 



10" 0.6 17.6 0.0012 0.27 0.03 

10i° 7 16.7 0.011 0.22 0.03 

10" 89 20.7 0.11 0.21 0.02 

4-10" 374 19.7 0.35 0.15 0.02 

6-10" 585 15.8 0.51 0.16 0.014 

10i2 1013 21.2 0.40 0.06 0.001 



different stellar population from which novae originate. 
This suggestion is supported both by observational and 
theoretical grounds. From one side Duerbeck (1990) and 
Delia Valle et al. (1992, 1994) have demonstrated, on the 
basis of galactic and extragalactic nova observations, the 



existence of two populations of novae: fast and bright no- 
vac belonging to 'disk' stellar population and slow and 
faint novae which originate from a 'bulge' stellar popula- 
tion. In particular, the latter authors (see also Delia Valle 
& Livio 1998) suggested that bulge novae could originate 
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from relatively light WDs, likely in the range of masses of 
{Mwd) < O.QMq while novae in the disk arise from mas- 
sive WDs {{Mwd) > l-^o)- From the other side, Truran 
(1990, see also Ritter et al. 1991) has found that the mass 
of the WD and the recurrence time between the outbursts 
arc inversely proportional. In order to match the 'theoret- 
ical' with 'empirical' rates, one needs to lower either a or 
n or both. For example for M87 the match between pre- 
dicted and observed rates can be achieved by increasing 
the recurrence time between two consecutive outbursts, to 

~ 3 — 1 • 10^ yr. This is about 10-3 times larger than 
assumed for novae in the 'disk' of the Galaxy. In this way, 
both the predictions of models I and II, for a luminous 
mass in the range 4-6 -W^^Mq are in very good agree- 
ment with the preliminary new estimate from Shara & 
Zurek (2002) of the nova rate for M87. In fact, M87 can 
be modeled both as an elliptical galaxy with luminous 
mass ~ 4 • 10"^"'^ Mq and Reff ^ 6 kpc or with luminous 
mass ~ 6 • 10^^ Mq and Reff ~ 7 kpc (Cohen & Ryzhov 
1997). The predictions of these two models do not differ 
much and both of them can well represent M87. Another 
possibility is to adopt a fraction of binary systems giving 
rise to novae smaller than 0.01. Unfortunately, if from one 
side we can justify this value for the Milky Way because 
this a gives both the observed present time Galactic rate 
and the correct yields of Lithium and CNO isotopes (see 
Romano & Matteucci, 2003), on the other side, we com- 
pletely lack this information for ellipticals. 

Finally, a comparison between tables 6 and 7 with 
Table 2 shows a very good agreement between the ob- 
served and predicted type la SN rates both in units of 
SNub and SNuk- 

In figures 2 and 3 we show the predicted star formation, 
nova and SNIa rates as functions of the cosmic time for 
ellipticals of different luminous masses and for the two 
IMFs. 

Figure 4 shows the nova rate in the ellipticals as a 
function of the mass of the galaxy. From this figure one 
can infer the expected nova rate as a function of luminous 
galactic mass. 

2.3. Magellanic Irregulars 

As an example of Magellanic irregular galaxies and late- 
type spiral we computed a model to fit the LMC. In par- 
ticular, we assumed an history of star formation with two 
main bursts occurring during the first 3 Gyr of the galactic 
life and during the last 2 Gyr, respectively, coupled with 
a lower star formation rate in between the two, as shown 
in Figure 5. This star formation regime for the LMC has 
been suggested by Calura et al. (2003) and it can repro- 
duce the observed [0/Fe] versus [Fe/H]. The star forma- 
tion rate is like in Eq. (5) with an efficiency of star forma- 
tion u = O.lGyr"^ in each starburst and oiv = O.OlGyr"^ 
in the interburst period. The IMF is the Salpeter one and 
the parameters a = 0.01 and A = 0.18. The history of star 
formation has been chosen to reproduce [0/Fe] vs. [Fe/H] 




5 10 15 

t (Gyr) 



Fig. 1. Nova rate and SNela rate predictions for the "two- 
infall" model for the Galaxy. Note the oscillatory be- 
haviour of the star formation rate during the last 2 or 
3 Gyr 

in the LMC (see Figure 5 and Calura et al. 2003). We con- 
sidered a luminous mass of ~ IO^^Mq which is suitable for 
the LMC (e.g. RusseU & Dopita, 1992). We adopted the 
same values of a and n as for the Milky Way and found 
that the predicted present time nova rate is in excellent 
agreement with the observed one and is ^ 1.7 nova yr^^. 
The predicted B and K luminosities are Lb = '2 ■ W^Lq 
and Lk = 0.27 • lO^^L©, respectively. As a consequence, 
we predict a nova rate in units of LNSRk of ~ 6, in agree- 
ment with the data of Table 1. The predicted type la SN 
rate is ~ 0.31SNu, also in very good agreement with the 
data of Table 2. 

3. Discussion and Conclusions 

In this paper we have computed the Galactic and extra- 
galactic nova rates, in particular the nova rates expected 
in elliptical galaxies and in Magellanic irregulars, typified 
by the LMC. To do that we have adopted detailed mod- 
els of chemical and photometric evolution of the Milky 
Way, ellipticals and the LMC. These models can repro- 
duce the majority of the observational constraints, includ- 
ing abundances, abundance ratios, luminosities and col- 
ors. We have computed also the type la SN rates since 
the novae and type la SNe are originating by CV-like bi- 
nary systems. The nova rate is calculated in a simple way 
but still depends on several unknown parameters such as 
the fraction of binary systems which are nova progenitors 
(parameter a) and the total number of outbursts occur- 
ring during a nova lifetime (parameter n). For the sake of 
simplicity, we have assumed these two parameters to be 
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Fig. 2. Ellipticals: star formation rate, SNela rate and no- 
vae rate predictions for model I. Dotted, lO^M©; short- 
dashed, 101"Mq; long- dashed, lO^M©; dotted - short- 
dashed, 10^2^0. 



Fig. 3. Ellipticals: star formation rate, SNela rate and no- 
vae rate predictions for model II. Dotted, lO^M©; short- 
dashed, 101°Mq; long- dashed, lO^^Mg; dotted - short- 
dashed, W^Mq. 



the same in all galaxies whereas the star formation his- 
tories are different and tuned to reproduce the available 
observational constraints. For the Milky Way this particu- 
lar choice of parameters, besides reproducing the observed 
nova rate, allows to predict "^Li, ^^N, ^^C and ^^O, in good 
agreement with the observations (Romano & Matteucci, 
2003). On the other hand, for external galaxies we do not 
have other constraints related to novae except the nova 
rates. Therefore, this approach docs not allow us to put 
strong constraints upon a and n singularly but rather on 
the product a ■ n. Concerning the delay with which no- 
vae appear since the beginning of star formation, we find 
that an average delay of 2 Gyr can be assumed for all the 
studied galaxies. 

Our main conclusions can be summarized as follows: 



— We predict that the nova rates in ellipticals reach a 
maximum after ~ 10^ years since the beginning of 
star formation and then they decrease after that. The 
same behaviour is shared by type la SNe although they 
reach a maximum after only a few 10^ years (see also 
Matteucci & Recchi, 2001). The delayed maximum in 
the nova rate relative to the SNIa rate is due to the 
longer time delay assumed for novae to appear relative 
to SNela (see Canal et al. 1996). 

— Our results are suggesting that a model for ellipticals 
which assumes a strong burst of star formation occur- 
ring at early cosmic times can reproduce the majority 
of the observational constraints in ellipticals. The same 
model predicts that the nova rates are increasing with 
the galactic mass of ellipticals and for a galaxy like 



9 10 11 12 

Fig. 4. Ellipticals: nova rate as a function of M;„m of the 
galaxy: model I (dotted), model II (short-dashed). 

M87 we predict a present time nova rate of 300-1000 
nova yr~^, when the values of a and n are the same 
as for the Galaxy. A good agreement with the obser- 
vations (100-300 nova yr~^) is obtained if we assume 
an average recurrence time between two consecutive 
nova outbursts of the order of 1 — 3 • 10^ yr, which is a 



8 



Prancesca Matteucci et al.: Modelling the nova rate in galaxies 





^ — : 




- 



-1.5 -1 

iFe/H] 



P 10 

5 

1 

k 2 

o 

^ 

03 

I -2 

to 



s° 0.5 



t (Gyr) 



t (Gyr) 



t (Gyr) 



Fig. 5. Irregulars: star formation rate, SNcIa rate, nova 
rate and [0/Fe] vs. [Fe/H] predictions for the LMC. The 
data are from Hill et al. 2000. 



factor of 10 larger than the typical recurrence time for 
novae originating in the disk of the Milky Way. This 
fact is a consequence of the existence of two popula- 
tions of novae, one associated with disk/spiral arms 
stellar population and the other with bulge/thick disk 
stellar population (Delia Valle et al., 1994; Yungelson 
et al. 1997). 

— The nova rate in ellipticals is rather independent from 
the adopted IMF, but it depends on the assumption 
that the cfhcicncy of star formation is an increasing 
function of the luminous galactic mass (inverse wind 
model of Matteucci 1994). We do not predict any trend 
for the luminosity specific nova rate as a function of 
the galactic mass for these galaxies. 

— We predict also a type la SN rate for ellipticals in 
very good agreement with the observed estimates and 
a slight decreasing trend for the SNIa rate with galactic 
blue and red luminosity. 

— Finally, we confirm the existence of a trend for the 
luminosity specific nova rate, as derived from the ob- 
servations (Delia Valle et al. 1994), passing from a 
Magellanic irregular (spiral) like the LMC, ~ 5, to a 
Sb spiral galaxy like the Milky Way ~ 2.1. For a giant 
elliptical like M87 we obtain a large range of uncer- 
tainty between 2.5 and 7.5 novae/yr per unit of 10^°K 
luminosity. The lower limit of the rate is very similar to 
that obtained from ground based nova survey whereas 
the upper limit is consistent with nova rates obtained 
from HST observations (Shara & Zurek 2002). If true, 
this could indicate that the studies of novae in extra- 
galactic systems from ground-based telescopes, can be 
largely affected by observational bias. As an alterna- 



tive, we note that Livio et al. (2002) have suggested 
that M87 could be a fast producer of novae, because 

of the presence of the jet which can enhance the accre- 
tion onto white dwarfs and therefore increase the nova 
rate. 
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